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£\\ , Abstract 

In a series of papers with Kossakowski, the first author has examined properties of densi- 
ties for which the positive partial transposition (PPT) property can be readily checked. These 
Oh| densities were also investigated from a different perspective by Baumgartner, Hiesmayr and 

Narnhofer. In this paper we show how the support of such densities can be expressed in 
terms of lines in a finite geometry and how that same structure lends itself to checking the 
necessary PPT condition and to a novel sufficient condition for separability. 
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1 Introduction 

7— I ■ 

£■ — ' Interest in quantum information theory has dramatically increased as the effectiveness of us- 
ing quantum entanglement as a resource for storing and manipulating information has become 
more apparent. Early ideas for applications include Shor's and Grover's algorithms (see [1] for 
example) as well as quantum key distribution [2], and those insights have motivated the surge 
■ in theoretical and experimental work during the last fifteen years. Since much of quantum in- 
formation theory is modelled in the context of finite dimensional composite systems, a subject 
of particular interest is determining the presence or absence of entanglement in terms of the 

. £h ' structure of the density matrix modelling the system. 

^ . If one is dealing with two d— level particles, the bipartite context, one can model the state of 

the composite system as a positive semi-definite trace one matrix p on the tensor product product 
space Cd®Cdi where Cd denotes a d-dimensional Hilbert space over the complex numbers. The 
system is said to be separable if it is in the convex hull of densities of the form t\ <S)T2, where each 
Tfc is a d-dimensional density on its respective space. An early observation by Peres [3] was that 
such separable states have the positive partial transpose PPT property: they remain separable 
under the partial transpose defined, for example, by r\ <g> T2 — > T\ <g) r|, where the superscript 
denotes the transpose operation. 

Since separable systems can be modelled classically, the two particles are not entangled. 
Thus the PPT condition is a necessary condition for separability, and it has been shown that it 
is also sufficient in the tensor dimensions 2® 2 and 2® 3 [4j. In higher dimensions, the condition 
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is not sufficient, as illustrated by a number of examples such as those in [5] and [6] among 
others. Nonetheless it has proved to be a surprisingly useful criterion, and there are a variety 
of examples of densities illustrating that fact. In [J] (see also [8]) the authors observed that 
many of these examples have a common property that facilitates checking the PPT property. 
That same structure was also defined in [9] , [TO] , and [11] , where the goal was to investigate the 
geometry of a particular subset of densities using "discrete phase space" as a tool. That discrete 
structure has been used extensively in other contexts - [12], [13], and [13] are just a few of many 
examples - and the first part of this paper relates the structure of this class of densities to the 
finite geometry of phase space and to checking the PPT condition. 

In the second part of the paper we show how the structure of these densities permits the 
development of a useful sufficient condition for separability. This condition is an offspring of 
a sufficient condition for separability that appears in [6], for example, and that serves as a 
counterpoint to the necessary PPT condition, since there are separable densities that don't 
satisfy it. The third part of the paper consists of a number of examples, including some of those 
in [3 [8] , [9] and [10] , and illustrates the use of this structural criterion 

We should note in passing that the investigation of separability and entanglement is an active 
area of research and there and too many papers to cite. Two survey papers on separability and on 
entanglement are [15J and [16J. The interested reader is referred to those surveys for an overview 
of the subject and also to [6] for the derivation of the sufficient condition that motivates our 
main result. 

2 Circulant densities and the PPT condition 

Motivated by the examples in [7], we establish notation that expresses the support of circulant 
densities in terms of lines in the context of finite geometry. The densities p to be considered are 
d 2 x d 2 positive semidefinite, trace one matrices with a particular pattern of support, and we 
begin by distinguishing between the positions where non-zero entries can appear and the entries 
themselves. Let M be a d 2 x d 2 matrix with entries equal to zero or one and let support(M) 
denote the set of positions of M with entry equal to one. 

Let {B (j, k) : < j, k < d} denote the constituent d x d submatrices, so that when d = 3, 
for example, 

/ £(0,0) 5(0,1) B(0,2) \ 
M= 5(1,0) 5(1,1) B(l,2) . (2.1) 
V 5(2,0) 5(2,1) 5(2,2) / 

For M to provide the locations for the non-zero entries of a density, we obviously need symmetry, 
M = M*, and that is equivalent to 5* (j, k) = B (k,j) for all j, k. Entries of M can be indexed 
as M(r, s), < r, s < d 2 or in tensor product notation as Mj 1 j 2 ,k 1 ^ 2 , < ji,j2> &l> &2 < d. The 
relationship between the two is 

(r, s) «-> (t#i + 32, dh + k 2 ) . (2.2) 

In particular, we routinely use the fact that such an entry corresponds to the (^'2,^2) entry of 
B(hM)- 
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As the authors noted implicitly or explicitly in [7J and |10j . the support of M can be inter- 
preted as lines in a two dimensional module; that is, they are lines in 



V 2 (d) = {(x,y):x,yeZ d } , (2.3) 

where Z d denotes the integers modulo d. (If d is prime, V 2 (d) is a vector space.) A typical 
example is the class of lines of the form 

L p = {(x, y) : y = x + p, x G Z d } , (2.4) 

where the addition is mod d, and one connects such lines with the matrix examples in [7J by 
orienting the y-axis down. For example, when d = 3 and p = 1, L\ can be represented by the 
ones in 

1 

1 
10 

We can also include "vertical" and "horizontal" sets of lines using L p to denote {(x,y) : x = p} 
or {(x,y) :y =p}. 

By analogy with the Euclidean plane, two lines are called parallel if they do not intersect. 
For two lines in the same class, L p n L q = <p if p ^ q, and it is easy to confirm the following 
result. 

Lemma 1 {L p : < p < d} is a partition ofV2(d) by a set of d mutually parallel lines. 

Note that we haven't defined lines using the more general formula ax + by + c = 0. This is 
because we are not assuming that d is a prime, and thus we cannot use the properties of an 
algebraic field which guarantee multiplicative inverses of non-zero elements. 

With this notation in hand, we can describe the pattern of non-zero entries in a generalized 
circulant density as developed in [7]. Let p denote a permutation of Z^, with the proviso that 
p (0) = 0. Define the non-zero entries of B p (j, k) by the line 

L p (j,k) ={(x+p(j),x+p(k)) :xeZ d } , (2.5) 

and note that B p (j, k) = B p (k,j) (When p is the identity permutation, we suppress the sub- 
script.) Then M p itself is defined in terms of Dirac notation and the L p (j, k) as 



M p = J2 

j,k 



d-1 



Y^\j) (k\®\p (j) + x) {p(k)+x\ 



x=0 



£/ p (x) , (2.6) 



where I p (x) = Ylj,k \ j) W ® \p U) + x ) (p (&) + x \ ■ 

As an example, let d = 3 and let p denote the identity permutation. Using x^ to denote an 
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entry of one when x = k and dots for entries of zero, the resulting M matrix is 

/ X ■ . ■ X Q . . ■ Xq \ 

X\ . . X\ X\ . 

X 2 X 2 ■ . x 2 ■ 

x 2 x 2 ■ . X 2 . 

X . . . X . . . X Q 

X\ . . X\ X\ . 

X\ . . X\ X\ 

x 2 x 2 . . . x 2 . 

\ X . . . XQ . . ■ XQ J 

Notice in particular the arrangement of the non-zero entries in the different 3x3 blocks: in each 
B(j, k) there is precisely one point of support associated with each x^. 

In confirming the properties of a putative density p, it is easy to check the trace one and 
the Hermitian conditions. What is harder is checking that p is positive semidefinite. One of 
the points of the present discussion is that if p has support in support (M p ), then that diffi- 
culty is reduced by the following observation. Representing a d 2 vector v in Dirac notation as 
Y^j Vjk \j) \k), it is easy to confirm that the entries of p are partitioned into disjoint sets indexed 
by x in the expression: 

H p \ v ) = ^^jix+pimPjix+pu^^ix+pik^Vkix+pik)) ■ (2-7) 

x j,k 

It is also easy to check that the components of v are partitioned into classes indexed by x, and 
thus one can check that p is positive semidefinite by checking that p restricted to each of the 
I p (x) is positive semidefinite. This means one is checking d d x d matrices rather than one 
d 2 x d 2 matrix. 

That same feature simplifies checking the PPT property, as was established in and |10j . 

Theorem 1 Let p denote a permutation of with p(0) = 0. Let M p be defined as in 112.6)) . 
Then M p is the sum of the disjoint matrices {L p (x) : x £ Z^} and for each (j, k) and each x there 
is exactly one index in support(B p (j, k)(lL p (x)). Lf p is a density with support(p) C support(M p ), 
then p is PPT if and only if p restricted to L_ p {y) is positive semidefinite for each y. 

Proof: The third sentence follows from the definitions. For the last assertion, let p T denote 
the partial transposition of p so that 

pT = \ x + P ( k )) Pj(x+p(j))Mx+p(k)) ( k \ {x+p{j)\ ■ 

j,k x 

Substitute x = y — p(j) — p (k) and reorder the summations: 

P T = /Z/~2\j) \V ~ PU)) Pj(y-p(k)),k(y-p(j)) ( k \ ( x ~P( k )\ ■ 
y j,k 

We see that the entries of p are again partitioned into disjoint sets indexed by y, and the earlier 
discussion applies, using the permutation —p, and completing the proof. □ 
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The 9x9 example above illustrates the idea. In tensor indexing, one of the three 3x3 
matrices of a density p and one of the 3x3 matrices of p T that one would check for positive 
semidefiteness are 

(/>00,00 POO, 11 P00,22 \ / P00,00 />02,10 P01.20 

Pll,00 Pll,ll Pll,22 I P T '■ I Pl0,02 Pl2,12 Pl2,20 

P22,00 P22A1 P22,22 / \ P20.01 P20,12 P21.21 

A useful necessary and sufficient condition to verify p is positive semidefinite is that the principal 
leading minors are non-negative. We leave it to the reader to write out the other submatrices 
for p and its partial transpose. 

In the preceding discussion we used addition of indices in the context of When d is the 
power of a prime, we have another alternative, and that is to use the Galois field GF{p n ). For 
example when d = 4 = 2 2 we could have defined the four parallel lines (right-side up) as 

/ A + l A 1 \ 

A A+l 1 

1 A+l A 

\ 1 A A + l / 

where GF{2 2 ) is the set {0, 1, A, A + 1} and the table above actually defines addition. (Multipli- 
cation uses A 2 = A + 1 and the usual properties of a field.). Since our only algebraic requirement 
so far is that addition be commutative, we could have used this different set of scalars in the 
prime power case. With the appropriate notation to index rows and columns, this gives an 
example of a d = 4 pattern that generalizes the circulant notation and that is not included in 
the appendix of [Tj- The proof of Theorem (2.1) then goes through unchanged. 




3 Generalized spin matrices and separability 

The PPT condition is a necessary condition for separability and is expressed in terms of the 
components of a density in the computational basis. Using a generalization of the Pauli spin 
matrices, one can obtain a sufficient condition for separability of a density in terms of the 
coefficients of a particular orthogonal family of unitary matrices. An accessible reference is 
|18j , and since the geometry of the support of the densities considered in this paper make them 
amenable to this approach, we summarize the salient points here. In [9] and [TO] this family of 
matrices is used to define entangled projections, and the papers analyze the convex hull of those 
projections. 

Let S = {Sj.k ■ < j, k < d} be defined by 

d-l 

S j)k = ^ rf m \m) (m + k\ , (3.1) 

m=0 

where r\ = e 2m / d . One can interpret these matrices as discrete Fourier transforms of the compu- 
tational basis matrices \m) (m + k\ , a feature that manifests itself in various applications. One 
always has Soo equal to the identity matrix. When d = 2, 5io = cr z , Sol = &x, and S\\ = ia y 
in terms of the usual Pauli spin matrices. That accounts for also calling them generalized Pauli 
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matrices, and for historical reasons they are also known as (discrete) Weyl-Heisenberg matrices 
or, possibly more accurately, as the discrete Weyl matrices [9]. 
Two useful properties of these matrices are: 



°j,k J u,v — '/ J j+u,k+vi 



\1 



n 



jk S _ 



j,~k 



o-l 
*j,k 



(3.2) 



Using the trace norm as an inner product, {A \B) = Tr (A*B), it follows that S is an orthogonal 
family of unitary matrices. As such, S is a basis for the space of d x d matrices, and a density 
can be written in this basis as 



1 

P= d 



j,k 



(3.3) 



where s ,o = 1 and s j:k = Tr S] k p 



. Prom earlier work we have the following result. 



Theorem 2 fWjj If p is a bipartite density and S u>v denotes the tensor product S UltVl (g) £, 
then a sufficient condition for p to be separable is that J2 U v \ i 
of S UtV in the spin representation of p. □ 



< 2, where s u v is the coefficient 



This is a relatively strong condition, although it is shown in [6] that the condition is sharp for 
certain families of densities with d = 2 n . The relevance here is that we can use the ideas behind 
the proof of this theorem and obtain sufficient conditions for separability for certain families of 
PPT circulant densities. 

To see the connection, let p denote the identity permutation. Then substituting m = u , 
n = m + t and v = m + k we can rewrite M as 



M 



E 



'rf-i 



\u) (v\ <g> \u + 1) (v + 1\ 

t=0 

|m) (m + k\ ® \n) (n + k\ 

k m n 
y^ •S'q,*; ® So,k ■ 



Thus the support of M equals the support of the sum of the tensor products £o,fc ® So,k- Since 
the support of Sj t k is the same as the support of 5o,fc, it is reasonable to assume that the spin 
matrix representation for a density with support (M) requires only tensor products of the form 
Sj lt k <S> Sj 2 ,k- That is indeed the case, and we skip the easy verification that the other spin 
coefficients equal zero. 

When the permutation p is not the identity, a similar analysis can be made that requires 
some additional notation, but no new concepts. We defer that discussion to the Appendix and 
stay with the assumption that p is the identity permutation. 

The next thing to notice is that the main diagonal of M is disjoint from any Sj u k ® &j2,k 
with k 0. Put another way, if po denotes the diagonal of a density p with support(M), then 



P~ PD 



<p S hk,32kSj 1 ,k &> Sj 2> k 

k^O 31,32 



(3.4) 
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Our strategy is to express the spin coefficients in terms of the computational basis entries of p 
and the spin matrices themselves in terms of projections. In certain cases we can then write p 
as a sum of tensor products of projections with non-negative coefficients - and that satisfies the 
definition of separability. To be effective we require that d be prime power and for computational 
simplicity we assume that d itself is prime. 

We provide the calculations in the Appendix and record here the basic structural result. 



Theorem 3 Suppose d is prime and support(p) is contained in support(M) as defined in A2.6\) . 
Then if po denotes the diagonal of p, there are rank one projections P a \ (m) such that 

d-l d-l 

P~PD= ^2 X] P ai ,i{mi) ® P a2j i(m 2 )C (a,m) , (3.5) 

ai,O2=0 mi,m2=0 

where each C (a, m) is real and d 2 C (a, m) equals 

E -(5')(ai+a 2 )-l;(mi+m2)-fc(oini+02ii2) ( c iR\ 
V K2J / J nxn2,(ni+fc)(n 2 +fe) ■ \ 6 - b ) 

4 First application of the structural result 



Equations (|3.5p and (j3.6h look a bit daunting, but in special cases they simplify quite nicely. 
Specifically, assume that p ni n 2 ,(n 1 +k)(n 2 +k)depends only on r = ri2 — n%: 

Pm{m+r),(ni+k)(m+r+k) = c r ■ (4-1) 

Then d 2 C (a, m) equals 

C r ^2 ??~^) (ai+a2)-k(m 1 +m, 2 +a 2 r) ^-ni(*(<Xl+a 2 )) ^ 2) 
r k^O ni 

= d<Ka 1 ,-a 2 )XAX) 7 r fe(mi+m2 ~ mi) 

= dS (a x , — Q2) Cr [dS (wi + TU2, ra L ) — 1] . (4-3) 

r 

4.1 A special case 

As an example, assume that c r = for r ^0: 

Pnin 2 ,(ni+fc)(n 2 +fc) = c5(ni,n 2 ) . 

Then 

P = Pd + ^ |n) \n) (n + k\{n + k\ , 

so the support of p is on the diagonal and on 1(0), in the notation above. (Recall that we have 
assumed d is prime.) We find that 

Q 

C (a, m) = -5 (ai + a 2 , 0) [dd (mi + m 2 , 0) - 1] . 
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Substituting this in (13, 5ft gives 

d-l d-l 

P = PD-^Y1 Yl P ai,l( m l)® P -ai,l( m 2) (4.4) 
ai=0 m\,m<2 
d-l d-l 

+C P ai ,l (mi) (8) P-oi.l (-mi) . 

ai=0 mi=0 

If c > 0, we can rewrite this using (16. 2ft to obtain 

d-l d-l 

P = PD - cl d <g> Id + c ^ ^ P ai ,% (mi) (8> -P-01,1 (~mi) . (4.5) 

ai=0 mi=0 

Thus the density /) will be separable provided the smallest entry in pp is at least c. 
Now suppose that c < 0. Then we can rework (|4.4p to obtain 

A> = PD-|c|(d-l)/ d ®/ d +|c|J^ P ai ,i (mi) ® P_ 0lj i (m 2 ) , (4.6) 

and /? is separable if the smallest entry in p£, is larger than |c| (d — 1). 

4.1.1 Isotropic state (example 2 from |8j) 

In this example c = X/d and the diagonal entries are either X/d + (1 — A) /d 2 or (1 — A) /d 2 . 
Then a sufficient condition for separability is A < 1/ (1 + d) , which also happens to be the PPT 
condition. 

4.1.2 Werner density (example 1 from [8j) 

Suppose one starts with a density of the form 

P = b h k ^ ^'1 + Y c J k 1-^) 0' fc l j 

so that the second expression on the right has only entries on the diagonal and the first term 
on the right includes all of the non-zero off-diagonal entries. For separability, it would suffice 
to prove either p or its partial transpose p T is separable, and it is convenient to work with the 
latter. Define 

1 ~P , P 
d+1 d-l ' 

The assumption in [8] is that for j ^ k, bji- = ac_, and bjj = x_ + x+. The diagonal entries 
Cjfc = x+. If x_ > 0, the sufficient condition trivially holds and the density is separable. If 
£_ < 0, the sufficient condition for separability from (|4.6p leads to p < 1/2, which again is the 
PPT condition. 



1 

X± = d 
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4.1.3 DiVincenzo et al example |17j (example 3 from [8]) 

The off-diagonal terms of p are ((c — b) /2) \jk) (kj\ in this example, while the diagonal terms 
are either (b + c) /2 or a = 1/d — (b + c) (d — 1) /2. (There is a misprint for the latter value in 
[8].) Again, the notation suggests working with p T to put it into the notational context of this 
section. If b < c, we find the sufficient conditions for separability are q 

< b and cd 2 + bd (d - 2) < 2 . 

If c < b, one obtains < c and bd{d—\) < 1. It is easy to check that these conditions are 
equivalent, so that 6 and c are non-negative in both cases and both of the other two conditions 
hold. As shown in [8], those were also the PPT conditions, so the separability analysis gives a 
stronger result than the PPT condition. 

4.1.4 Horodecki example (example 4 from |8j) 

In this example, d = 3 and once again n\ = ri2 in the notation above is a necessary condition 
for the support of p. c = 2/21, and the diagonal entries are either a/21 or (5 — a) /21. Then 
separability is guaranteed by 2 < a < 3 while the PPT condition allows 2 < a < 4. The theory 
used here does not resolve the question of separability when 3 < a < 4, and it is known that in 
this case the state is actually entangled. 

4.2 The general case 

We assume the c r 's are real and define s = ^ r c r . If we substitute into (|3.5j) . we can obtain 

p = p D - sl d <g) I d + ^2 c r ^2 X] Pa > 1 ( m ) ® jP ~ a > 1 ( ar ~ m ) • ( 4 ' 7 ) 

c am 

An easy case is when each of the c r 's is non- negative, since that gives 

s < min (p nk!nk ) , 

as a sufficient condition for separability. It is interesting to note that for d = 3, if we have 
\J^2k c k — m ™ (Pnk,nk) > then p satisfies the PPT condition and this is a weaker condition than 
the sufficient condition for separability. 

As another special case, suppose that co < < c r for r / and also that s = 0. Then 
following the usual approach we can write p as a sum of separable projections with non-negative 
weights and a term pd — d |co| Id ® Id, giving the obvious sufficient condition. 

4.2.1 Baumgartner, Hiesmayr and Narnhofer: |9], |10] 

Define = Sj k (8 Id J2u \ u ) \ u ) aim ^ ne projection Pj^ = • The class of densities 

studied in [9] and [TO] are defined by 

P = Yl c i,kPj,k , (4.8) 
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where the constants are non-negative and sum to one. One finds that 



Pn(n+r) ,(n+k)(n+r+k) 



i >*■ (4-9) 



d 

3 



There is no dependence on n, and one calculates 

C (a, m) = -^5 (ox, —02) } y Cj >r [dS (air — j, mi + m 2 ) — 1] . (4.10) 

3 r 

In recent work on discrete Wigner functions, there is a connection between lines in phase 
space and projections - see for example [13] and [19] among many others. This motivated one 
of the examples in [9], and the general case is 

l/d , if r = sj +t 
, otherwise 

where s and t are fixed. The geometry is that one is defining a line in the index set, i.e. a discrete 
phase space. This includes horizontal lines, and vertical lines have to be treated separately. We 
then have C (a,m) equals 

■^6 (ai , -a 2 ) ^2 \ ( aiS -1), mi + rri2- ait) - 1] . 
3 

When a\s 7^ 1, there is exactly one value of j for which the (^-function equals 1 and thus 
C (a, m) = 0. If s = 0, this is true for all a and m, and p is diagonal and therefore separable. In 
the remaining cases, let a\ = a s = s , so that C (a, m) equals 



■^5 (ai,s x ) 5 (a 2 , 



-1 + (ts 1 ,mi + m 2 ) 



We find 



p = PD~ -pld ^Id + ^y. P °*A ( m ) ® p -a s ,i ( a st - m) 



Since the diagonal entries are p n (n+r),n(n+r) = \ ^L c j,r = ^> it follows that these densities are 
separable. From properties developed in the Appendix, one can also check that they are rank d 
projections, as noted in [9]. 

It is known that the normalized identity is contained in an open set of separable densities. 
We can observe that property in this context by considering parameterized segments of the form 

, , 1-i 
P W = ~cP~ Id2 + tp ' 

where p is defined by (|4.8p . We will show that if t < t 2 = 1/ (1 + d) then p (t) is separable. That 
same bound holds for Werner densities in the bipartite case, and the bound t n = 1/ (l + d™" 1 ) 
works for multipartite Werner densities [20]. The obvious conjecture is that t n works for the 
multipartite versions of (I4.8D . 



10 



The proof is an easy application of the structure result. Ignoring terms with non-negative 
coefficients, look at the remaining part of p : 

PD ~ ^2^2 XI Pa ^ ® P - a > 1 ( m2 "> ^2 c i' r = PD~ 2 I d? ■ 

a mi ,7?i2 j,r 

The usual requirement takes the form 

l-t t ^ t 

3 

so that for all r 

t < 1/ M +d-d^2cj )T 

This verifies the sufficiency of t < t% for separability. 

Finally, we illustrate the structural results for one other example from [9] where d = 3. a 
and (3 are non-negative parameters and in that notation 

l-a- (3 

P = 7. 1& + CtPl,0 + PP2,0 ■ 

y 

Putting this in the notation above, t = a + /3, ci,o = a/t and C2,o = P/t. Then the sufficient 
condition for separability ist = a+/3 < 4. Equation(50) in [9] provides a necessary and sufficient 
condition for p to be PPT, and it is easy to confirm the condition for separability satisfies that 
constraint. 

5 A class of circulant densities with product entries 



In the examples above, the entries of p in (14. 1|) did not depend on n. In this section they do, 
but in a very structured manner: 

Pn{n+r),(n+k)(n+r+k) = x { n i r ) x ( n + &i T) . 

This seems to be a novel set of examples that is amenable to analysis via the structure results. 
The key is that 

r ^-{ k 2 )(ai+a2)-k(mi+m2)—kn(ai+a2)—kra2 

can be written as 

^(2)( a i+ a 2)+n(m.i+m2+ra2)^-( n ^ ft )(ai+a2)-(n+fe)(mi+jTi2+ra2) 

This meshes well with x (n, r)x{n + k, r) to give d 2 C (a, m) + ^ n ^2 r \x (n, r)\ 2 equal to 

2 



x ( nj r ) ^(2)( a i+ a 2)+n{mi+m2+r. 



C12> 
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Define 

A(r, b,t) 

and 



d 

n 



Q(r,b,t)= Pa i ,i{m 1 )®P a2> i{m 2 ) 

fe=ai+a 2 t=mi+m 2 +ra 2 

so that up to normalization the Q's are separable. Then the separable representation is 
P = PD-^\x{n,r)\ 2 Id®Id + ^A{rAt)Q{rAt) 

n,T <Ji,<J 2 r,b,t 

and it follows that a sufficient condition for separability is 

min {pmn.mn) > V" \x (n, r) | 2 . 

n,r 

Denoting one of these densities as p x we can define the line segment 

/ N 1 - * r 

Px (t) = —^2~ I d? + tPx 

and obtain a sufficient condition for separability 

(n, r)\ 2 - min (p mn , mn ) < 1 • 



i ( 1 + ( ^ \ x ( T 



rn,n 



We omit the details. 



6 Appendix: technical details of the proof of Theorem 3.2 

To obtain the representation above of p — po, we first need to relate projections to the spin 
matrices, and here is where we restrict the discussion by assuming that d is a prime. Fix an 
index (j, k) and define 

m 

Using the properties above, it is straightforward to check that 

P hk {r)P hk {s) = 5{r )S )P hk {r) , 

and that the family {Pj k ( r ) : < r < d} is an orthogonal family of rank one projections. (As 
an aside, for the reader interested in mutually unbiased bases, this is one way to construct them 
when the dimension of the space is a prime power. See [14].) Note that one could interpret (|6.ip 
as another discrete Fourier transform of matrices, so that one should be able to compute the 
spin matrices from the associated projections. In fact the formula is 

(v r Sa,iY = £ r mt Pa,x (m + r) , (6.2) 
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and this is nothing more than the spectral decomposition of ijf S a ,if- When t = the left side 
is interpreted as the identity. 

Now using (|3,2|) it is easy to confirm that (S a: i) k = rf\v) S a ^^ provided we interpret the 
binomial coefficient as if k = or k = 1. Since d is prime, for k ^ and given j we can define 
a = jk^ 1 so that 

% fe = rr<*)^rT mfc Pa,i(m) , (6.3) 

m 

and that completes the first step. 

For the second step we use the structure of the support of p to write 



m^v 



d-1 



Pm(m+t),v(v+t) \ m ) H ® |m + t) (v + t\ 



t=0 



= Yl Yl Pnm 2 ,{ni+k)(n 2 +k) («1 + H <g> |n 2 ) (n 2 + fc| . 

Comparing this with (|3.4p we conclude that the constants Pnin2,(ni+fc)(n 2 +fc) are transforms of 
the Sj 1 kj 2 k j 

A y nin 2 ,(ni+fc)(n 2 +A:) _ ^2 ' ' °Jlk,j2k ; 

ii -J2 

and thus that 

S JifcJ 2 fc = S S r ?~ niJl ~" 252 /'n 1 n 2 ,(n 1 +fc)(n 2 +fe) • (6.4) 
ni r» 2 

Substituting (|6.4p and (|6.2p into (|3.4p gives Theorem 3.2, up to the confirmation that C(a, m) 
is real. That last fact follows readily enough by taking the complex conjugate of C(a,m), 
substituting —k for k and then making some additional notational changes to obtain C(a,m). 

Finally, we show that assuming the permutation p is not the identity only complicates the 
notation without generalizing the discussion. Suppose then that the permutation p defining M p 
is not the identity. Let a denote the inverse permutation p~ l . Make the substitutions p (u) = m, 
n = m + t, and p(v) = m + k in the expression for M p to obtain 

M p = ^2^2^2\v){v\®\p(u) + t)ip(v) + t\ 

t U V 

= Y2 Y2 1* 7 ( a ( m + ^)i Y2 ( n ^ 

km n 

= ■S'o.fc (o" ) 'S'o.fc 



where Sj^ (a) = if m \a (m)) (a (m + k)\ . Introducing the permutation a in the definition 
of the spin matrices is equivalent to a relabeling of the computational basis in the first space 
and the class S (a) should have the same properties as S. It is easy to check that is indeed the 
case, and thus the only difference that one has in (13. 5p is that the projections P nii i(mi) should 
have a a added to the notation. Thus, for all practical purposes in studying classes of examples, 
we can simply assume that p is the identity. 
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